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SUMMARY 


RosultB of computer-code time-dependent aolutions of the two-dimensional 
compressible Navier-Stokes equations and the results of Independent experi- 
ments are compared to verify the Mach-number range for instabilities in the 
transonic flow field about a 14%-thick biconvex airfoil at an angle of attack 
of 0° and a Reynolds numbf^ of 7x10*’. The computer code predictions were made 
at Ames Research Center, <d the experiments were conducted in a transonic, 
slotted-wall wind tunnel &*• the Royal Aircraft Establishment, Bedford, England. 
The computer code included an algebraic eddy-viscosity turbulence model devel- 
oped for steady flows, and all computations were made using free-f light bound- 
ary conditions. All of the features documented experimentally for both steady 
and unsteady flows were predicted qualitatively; even with the above simplifi- 
cations, the predictions were, on the whole, in good quantitative agreement 
with experiment. In particular, predicted time histories of shock-wave posi- 
tion, surface pressures, lift, and pitching moment were found to be in very 
good agreement with experiment for an unsteady flow. Depending upon the free- 
stream Mach number for steady flows, the surface pressure downstream of the 
shock wave or the shock-wave location was not well predicted. 


NOMENCLATURE 


C„ lift coefficient 

pitching-moment coefficient about the airfoil quarter chord 
C,, pressure coefficient, (p - p„)/q„ 

r 

r airfoil chord 
M Mach number 

p pressure 

q dynamic pressure 

Re^ Reynolds number based on airfoil chord and free-stream conditions 

T temperature 



t 


time 


u velocity component in the x direction 

Xty axes parallel and normal, respectively, to airfoil chord with origin at 
airfoil leadin;^ edge 

u circular frequency 

Subscripts: 

«> free-strearo conditions 

rtns root-mean-square value 

t total conditions 

Superscript; 

* sonic conditions 


INTRODUCTION 


For a manlier of years Ames Research Center has been conducting companion 
programs of basic research in experimental and computational fluid dynamics 
for testing and guiding the development of turbulence modeling within regions 
of separated flows (ref. 1). One of these programs has contributed signifi- 
cantly to a renewed Interest in hlgh-Reynolds-number transonic flows, specif- 
ically the experimental documentation and computation of transonic flows about 
an 18%-thlck circular arc airfoil (refs. 2-9). For this airfoil at an angle 
of attack of O'* , three distinctly different types of supercritical flows were 
documented experimental ly tor the first time in reference 3. Similar documen- 
tation at angles of attack was subsequently reported In reference 8. For a 
constant angle of attack and Reynolds number, and at progressively higher 
froe-streara Mach numbers, the flow field about the airfoil can be (1) steady, 
with tral ling-edge boundary-layer separation, (2) unsteady, with aerodynami- 
eally self-excited periodic oscillations In shock-wave location and intensity 
and in the location of bound iry-layer separation, or (3) can be steady with 
shock-induced separation. It was demonstrated in reference 6 that two- 
dimensional computer codes for solutions of the Reynolds-averaged, compres- 
sible Navier-Stokes equations are capable of reproducing the time-dependent 
features of unsteady, transonic, turbulent separated flows, even if the code 
employs a simple algebraic eddy-viscosity turbulence model developed for 
steady flows. 

An independent experimental study of the effect of thickness-to-chord 
ratio on the transonic Mach-number range for shock-induced periodic flows 
about biconvex airfoils at an angle of attack of 0° and at Reynolds numbers 
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up to exlO*^ has bean made uy Mabey at the Royal Aircraft Eatabliahment . 
Bedford, England (ref. 10). Baaed on the reaulte of reference 10, the RAE 
decided to extend the study to Reynolds numbers up to 7x10 . For this pur- 
pose a 14%-thlck biconvex airfoil was selected to verify the expected Mach- 
manber range for transonic- flow Instabilities and to measure the chordwlt.’ 
distribution of oscillatory pressures. Furthermore, the results could provide 
additional experimental data against which to test the feasibility of the Ames 
Navler-Stokes-type computer code, used in reference 6, for pic-dictlng both 
steady and unsteady transonic, separated turbulent flows. Accordingly, NASA 
and RAE entered into an agreement in which Ames Research Center would compute 
three Navier-Stokes solutions to complement the experimental results given in 
reference 11. Using the results of reference 10 as a guide, cne solution 
would be at a Mach number to verify an expected oscillatory flow and one solu- 
tion each at lower and higher Mach nvanbers to verify expected steady flows. 

This paper briefly describes the experimental and computational 
approaches to this joint study. Similarities and differences between the 
experiments and computations are noted, and computed and experimental rosuJts 
are compared to illustrate the current predictive capabilities of a two- 
dimensional computer program coded for the Reynolds-averaged, compressible 
Navier-Scokes equations and an algebraic eddy-viscosity turbulence model 
developed for steady flows. 


AFPRC '’Jl 


Selection of the NASA/RAE-collaboratlon-agreement test conditions for 
comparisons of predicted and experimental results was based originally on the 
results of reference 10. Those results indicated that a 14%-thick biconvex 
airfoil at an angle of attack of 0° and at a Reynolds number of 6x10® exper- 
ienced unsteady periodic flows in the Mach-number range from 0.81 to 0.85. 

From the trends with Reynolds number indicated in reference 8, one would 
expect similar transonic-flow instabilities over the same small Mach-number 
range at a Reynolds number of 7x10®, the highest experimental Reynolds number 
attainable for the experiment (ref. 11). Accordingly, three sets of initial 
conditions were selected for the two-dimensional Navier-Stokes solutions. 

For a Reynolds number ROj, of 7x10®, free-stream Mach numbers M„ of 0.80 
and 0.86 were selected to verify expected steady flows; and an M„„ of 0.83 
was selected to verify an expected unsteady, oscillatory flow. The computa- 
tions for these three cases were initiated at Ames Research Center before the 
experiments were begun at the RAE. 

When the experiments were completed at RAE, the results indicated the 
desirability of changing M^o for the three computational cases because the 
experimental Moo range for periodic flows was found to be from 0.82 to 0.86. 
The free-stream conditions finally selected for comparison of predicted and 
experimental results are listed in table 1. Note that the computational free- 
stream conditions for the two steady-flow cases (M^o » 0.81 and 0.88) were 
changed to match the experiments. Computations at these Mj„'s were completed 
prior to receipt of experimental results. By the time the experiments were 
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TABLE 1.- FRBE-STREAM CONDITIONS 




RecXlO”* 

T,, *R 

Conputatlon 

0.81 

6.9 

530 

Exparlmant 

0.81 

6.9 

530 

Cooputafcion 

0.83 

7.0 

500 

Exparlnant 

0.85 

7.0 

532 

Computation 

0.88 

7.1 

534 

Experiment 

0.88 

7.1 

1 

534 


completed and detailed unsteady-flow data obtained at ■ 0.85, too much 
time had been Invested In the original calculations at Ma> ■ 0.83 to warrant 
a restart of this case, which exhibited an unsteaiy periodic flow. Other 
differences between the experiments and the computations are noted below. 


Experiment 

As stated in reference 11, an aluminum model of a 14%-thlck biconvex 
airfoil — chord 300 mm, semispan 600 ram — was mounted from the center of one 
side wall of a 640-mm-high by 910-iran-wide RAE wind tunnel with slotted upper 
and lower test-section walls. Although this arrangement is technically a 
three-dimensional test setup, the flow was determined to be "two-dimensional" 
at the center of the semispan model where both the static and dynamic pres- 
sures were measured simultaneously in real time at 15 chordwise stations on 
both upper and lower surfaces of the airfoil. To Insure a turbulent boundary 
layer over most of the airfoil at a Reynolds number of 7x10®, carborundum 
transition trips were used at 2% chord on each surface. The ventilated test 
section was unchoked during tests at all Mach numbers. 


Computations 

Predictions of the flow field were made using a slightly modified version 
of the computer program described in reference 6. That program utilizes an 
explicit finite-difference method (refs. 2, 12, 13) to solve the time- 
dependent, two-dimensional, Reynolds-averaged form of the compressible Navier- 
Stokes equations; it incorporates the efficient numerical solver described in 
reference 14. The computer code also includes the additional equations for 
the turbulence model described in reference 6. The slight modifications 
referred to above Include correction of the two subtle coding errors mentioned 
in reference 6, which earlier may have introduced numerical asymmetries. 

Also, the mesh size of the C-type computational domain was increased from 
78x35 to 138x35 to provide more mesh points on the airfoil in the streamwise 
direction. 
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Dctallad dtaerlptions of tha control voIum* aaah, and boundary condi- 
tlona ara glvan in rafaranca 15. For tha praaant atudy* a UX-thlek circular- 
arc airfoil waa Impulalvaly atartad fron rant at tlna laro at aach of tha 
ctmputatlonal fraa-atraaa eonditiona liatad in tabla 1. nia control voIum* 

-6 and 4-8 chorda in tha x diraction and t6 chorda in tha y diraction, waa 
dividad into a 138x35 vtah with 49 naah pointa on aach of tha uppar and lowar 
eirfoll aurfacaa. Tha flow-fiald davalof»Mnt within thia voluma waa followad 
in tine until, ncnninally* it attalnad a ataady atata. For tha praaant coda, 
airfoil, and fraa-atrawn eonditiona. thia required a tiaa equivalent to tha 
mean flow traveling about 8 chord langtha pact the airfoil (u»t/c > 8). Since 
tha diatributiona of flow velocity through tha alottad walla and the down- 
atream atatic preaaurea were not available fron tha axperinanta. fraa-atream 
eonditiona ware aaaumad at the far upatraan and tranavaraa boundariaa. and 
zero gradianta were aaaumad for all flow variablaa at tha dovnatraam boundary. 

To avoid numerical dlfflcultlaa uaaoclated with a aharp-edged airfoil and 
a C-type meah. all computationa ware made for a 14%-thlck circular arc airfoil 
modified for a 0.71% noae radiua. Thia foraahortena an unmodified airfoil 
approximately 2% and doea not affect reaulta computed at an angle of attack of 
0°. Consistent with the experiments of reference 11, turbulent boundary- layer 
flow was initiated over the airfoil at 2% chord. 


RESULTS AND DISCUSSION 


The computed results for the free-stream conditions listed in table 1 
predicted the same general features of the flow fields observed in the corre- 
sponding experiments reported in reference 11. At Ma, ■ 0,81 the flow field 
was steady with a strong shock wave and trailing-edge separation. At 

“ 0.83 the predicted flow field was unsteady, with periodic oscillations 
in shock-wave location and intensity and accompanying oscillations in the 
extent of boundary- layer separation. The same phenomena occurred on opposite 
surfaces of the airfoil 180® out of phase. Thia type of flow was observed 
experimentally at between 0.82 and 0,86. At ■ 0.88 the flow field 
was steady with a shock wave sufficiently strong to Induce separation at the 
base of the shock wave, which extended downstream beyond the airfoil trailing 
edge. The two steady flow fields, in actuality, were pseudosteady. The 
shock-wave location and intensity, the location of boundary separation, and 
the pressure distributions did experience small variations about mean values 
in both the computations and the experiments. 

Results of comparing details of the predicted and experimental data will 
be discussed separately for the steady and unsteady flows. To obtain root- 
mean-square values (rms) for the steady flows, predicted results for these 
cases were obtained from analyses of solutions computed over a time equivalent 
to 2 chord lengths of travel of the mean flow beyond that for which a solution 
is normally considered to have converged to a steady state. Because >f the 
symmetry of the mean quantities determined over this time interval, results 
for the steady flow fields are presented for only half of the airfoil. Pre- 
dicted results for the unsteady-flow case are from analyses of time histories 
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of tolutlons for four consacuclvc cyclca of flow f laid oacillatlon (approx- 
Imataly 27 chord langtha of travel of the nean flow). Certainly «ore cyulea 
are dealrable; however* four are conaiiered aufficient for purpoaea of >;he 
present compariaona. 


Steady Flowa 

Predicted and experimental ahock-wave locatlona and pressure dlatribu- 
tlona for 0.81 are ahown In figure 1. The Mach contoura In figure 1(a) 

Indicate that the predicted ahock wave la located leaa than 2 % downstream of 
the actual ahock wave. (Coaleacence of the contoura about the daahed aonlc 
line indicate the formation of a ahock wave the closer the contours* the 
stronger the shock wave). The same relative locatlona 'jre Indicated In fig- 
ure 1(b) by the predicted and experimental jumps In the mean pressure distri- 
bution near x/c “ 0.7. The mean pressure la predicted fairly well except 
over the last 20% of the airfoil. In the experiment* tralling-edge separation 
is reported ro occur at x/c ■ 0.87 (ref. 11)* while In the computations* 
separation occurs first at the base of the shock wave (x/c ■ 0.7) followed by 
reattachment and subsequent separation extending into the wake. The locations 
of reattachment and subsequent separation vary with time. An instantaneous 
example is indicated in figure 1(a) by the dotted line* which is the Mach 
contour defined for u ■ 0. 


-PREDICTED 
O EXPERIMENT, RiP. 11 



(a) Mean shock-wave location. 

Figure 1.- Predicted and experimental shock-wave locations and pressure 

distributions, Mo» ■ 0.81. 
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(;>) Mean pressure, 
(c) Rtns pressure. 
Figure 1.- Concluded. 


The distributi* a of the rms values of the broadband pressure fluctua- 
tions are compared in figure 1(c). The results are qualitatively similar in 
that an increase in pressure fluctuations occurs at and immediately downstream 
of the shock wave. The magnitude of the increase, however, is overpredicted. 
This overprediction is attributed to an unsteadiness in the computed location 
of the reattachment and subsequent trailing-edge separation (noted above) that 
was not observed in the experiments. 

The present computations, using free-f light upstream and transverse 
boundary conditions and zero-gradient-downstream boundary conditions, provide 
a good approximation to the flow field about this airfoil in the RAE slotted- 
wall test section at this M„o. This result is consistent with the fact that 
the sonic lino is computed to extend above the airfoil only to a y/c of 0.7 
(fig. 1(a)), about 34% below the slotted walls located at y/c •• 1.067. Dif- 
ferences in the predicted and experimental results noted in figure 1 are 
attributed to dif ficiencies in the turbulence model and a mismatch in tne 
downstream boundary conditions. 
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Predicted and cxparlnicntal ahock<-wavc locations and pressure distribu- 
tions for ■ 0.88 are shovn in figure 2. The predicted mean pressure is 
in excellent agreement with experiment ahead of the shock wave; however, as 
indicated by the Mach contours and pressure distributions in figures 2(d) and 
2(b), respectively, the shock wave is predicted to be 13% farther downstream 
than In the experiment. Prediction of the nearly correct mean pressure level 
downstream of the shock vave is fortuitous. In the experiment, the shock wave 
was nearly normal (see fig. 2(a)), with separation at the base of the shock 
extending downstream info the wake. The same type of separation pattern was 
predicted by the calculations (see, e.g., the dotted line in fig. 2(a)); how- 
ever, the separation emanates from the base of a calculated oblique shock wave. 

The above-noted differences may be due, in part, to inadequate turbulence 
modeling; however, they are believed to result primarily from the mismatch in 
boundary conditions between the computations and the experiment. For example, 
the sonic line was computed to extend 2 chord heights above and below the 
airfoil -- well beyond the location of the experimental slotted walls 
(y/c ■ 1.067). Although the location of the sonic line in the experiment is 
unknown, it is physically impossible that it extend much beyond the slotted 
wind-tunnel walls. Hence, there results a significant mismatch in the compu- 
tational and experimental tunnel-wall boundary conditions. In addition, with 
zero-gradient-downstream boundary conditions there is no assurance that the 
flow parameters computed at the x/c of the wind-tunnel exit plane would 
match the experimental parameters at that location, had they been measured. 


- PREDICTED 
J EXPERIMENT REP. 11 



(a) Mean shock-wave location. 

Figure 2.- Predicted and experimental shock-wave locations and pressure 

distributions, Moe, ■ 0.88. 
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(b) Mean pressure, 
(c) Rms pressure. 
Figure 2.- Concluded. 


Similar to the results at ■ 0.81, the rms pressure fluctuations 
Increase at and immediately downstream of the shock wave; however, at 
Mo, “ 0.88 the magnitude of the Increase is underpredicted except near the 
trailing edge (fig. 2(c)). Overprediction at the trailing edge is attributed 
to an intermittent separation and reattachment of the boundary layer computed 
in this region but not observed in the experiments. 


Unsteady Flows 

As noted earlier, both the predicted and experimental flow fields exhib- 
ited periodic oscillations in shock-wave location and Intensity, and in the 


9 



location of boundary- layer saparation on oppoalca surfacaa of the airfoil 160” 
out of phaia^ Although detailed conpariaone can be nade only between data 
available at different free-etream Mach numbere (aee table 1), it ia aignifi- 
eant to note that the reduced frequency of the flow-field oacillation wc/uh 
waa predicted to be 0.934; that reported in reference 11 at the aame Mach 
number (M«, • 0.83) w» 0.946. 

Predicted and experimental preaaure diatributiona for the unateady-f low 
caaea are compared in figure 3. In apite of the differencea in M^ in the 
two aeta of data, the experimental diatributiona are predicted quite well, 
particularly the dlatribution of rma preaaure fluctuatlona in figure 3(b). 
For thia caae it ia believed that Che effect on the computatlona of the mla- 
match in boundary condltlona ia far leaa Chan that at Mo, > 0.88. Thia la 
becauae examination of Mach-contour plota for a typical cycle of flow-field 
oacillation revealed ci.r.t the aonic line on each aide of the airfoil reached 
valuea of y/c Chat correspond to the location of the slotted tunnel walla 
during only 10% of Che cycle. Also, aubarquent comparlaons of time historieB 
of various aerodynamic quantities show good agreement between predicted and 
measured values. 

Time hlacoriea of shock-wave location during a typical cycle of oacilla- 
clon arc compared in flgsre 4(b). Shown in figures 4(a) and 4(c) are plots 
of computed Mach contours selected to Indicate the relative strength of the 
upper surface shock wave at its extreme positions during a cycle of motion. 

The "leaders" to figure 4(b) Indicate the times during the cycle for each 
contour plot. In order to minlmlre the effect of different free-stream Mach 
numbers, the shock positions (and all subsequent time histories) are compared 
on the basis of chords traveled by the mean flow (u^^t/c). The experimental 
results were inferred from time histories of surface-pressure distributions. 
The predicted shock positions were defined as the values of x/c at the mid- 
point of the jump in surface pressure coef flcien*; from 25 flow-field solutions 
used to define one cycle of oscillation. As an eyample, the data in fig- 
ure 2(b) indicate a shock wave on the surface locat^jd at x/c ■ 0.816 where 

Cp - -0.6. 


It was noted in reference 11. and is indicated in figure 4(b) that as 
the shock waves moved rearward they diminished in strength and at the most 
rearward locations were not visible in the high-speed shadowgraph movies. 
Although the Mach contours obtained during these portions of the cycle do 
indicate a definite sonic line similar to the dashed line shown above the air- 
foil in figure 4(a), it can be noted that the Mach-number gradient across the 
sonic line is much less (larger spacing in the streamwise direction) than when 
the shock wave is at the more forward chord stations (see. e.g., upper surface 
in figure 4(c)). In fact, values of the second derivative of the density 
(the quantity to which a shadowgraph system is sensitive), computed across the 
upper surface sonic line as far away from the surface as y/c *0.3 for 
figure 4(a) was found to decrease by a factor of 27 within the first 10% of a 
chord above the airfoil compared with that for figure 4(c). This result indi- 
cates that the intensity of the shock wave does greatly diminish away from the 
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airfoil aarfaca aa it movaa raarward in tha cMiputationB and probably 
vaniahaa* aa in eha cxparinant; howavar, tha unknown ralationahip ia how tha 
photeacnaitiva asulaion of high-apaad aovia fllai would raapond to tha magni- 
tuda of tha conputad variation in danaity. 
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(a) Mean pressure. 

(b) Rms pressure. 

Figure 3.- Predicted and experimental pressure distributions for unsteady flows. 
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CHORDS TRAVELED 


(a) Upper shock wave, rearward. 

(b) Shock-wave locations. 

(c) Upper shock wave, forward. 

Figure 4.- Predicted and experimental shock-wave locations during a typical 
cycle of oscillation and selected Mach-contour plots. 
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Time histories of the pressure at four selected chord stations are com- 
pared In figure 5 for two cycles of oscillation. The waveform, frequency, and 
magnitude of pressure variations are predicted quite well at the more rearward 
x/c locations (figs. 5(c) and 5(d)) where the rms pressures are also better 
predicted. See, for example, figure 3(b) at x/c *0.7 and 0.8. The most 
forward movement of the shock wave shown In figure 4(b) is to x/c - 0.65. 

In the experiment, weak pressure fluctuations propagate as far forward as 
x/c • 0.55. as shown In figure 5(a) and 5(b). In reference 11 these fluctua- 
tions are attributed to weak compression waves ahead of the main shock wave, 
which apparently are not present In the computations. 

Even though the magnitude of th<>- instantaneous surface pressures are 
overpredicted and underpredicted at different chord locations, the chordwlse 
integration of instantaneous pressures generally overpredicts the magnitudes 
of the unsteady lift and pitching moment about the quarter chord. These pre- 
dicted and experimental time histories are compared in figure 6 where xt can 
be seen that the reduced frequency of oscillation is predicted quite well. 

PREDICTED, - 0.83 

EXPERIMEMT, REF. 11 



x/c • 0.70 x/c « 0.80 

CHORDS TRAVELED CHORDS TRAVELED 

Figure 5.- Predicted and experimental surface-pressure time histories. 
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PREDICTED, M,, - 1S3 

EXPERIMENT, M„- 0.86. RiP.11 



0 8 16 24 32 40 

CHORDS TRAVELED 


(a) Pitching moment. 

(b) Lift. 

Figure 6.- Predicted and experimental pitching-moment and lift time histories. 


CONCLUDING REMARKS 


A computer code for time-dependent solutions of the tv.'o-dimensional , 
Reynolds-averaged form of the compressible Navier-Stokes equations has been 
used at Ames Research Center to predict and compare results with experiments 
made at the Royal Aircraft Establishment, Bedford, England, to verify the 
Mach number range for instabilities in the transonic flow field about a 


14 



I4%”thick biconvex airfoil at an angle of attack of 0“ and a Reynolds number 
of 7x10®. The experiments were conducted in a transonic wind tunnel with 
slotted upper and lower walls. The computer code included a simple algebraic 
eddy-viscosity turbulence model developed for steady flows, and all computa- 
tions were made using free-fllght boundary conditions. All of the features 
documented experimentally for both steady and unsteady flows were predicted 
qualitatively, and in many instances the predictions were in good quantitative 
agreement with experiment, even with the above simplifications. 

At a free-stream Mach number of 0.81, where the characteristic features 
of the steady flow are a reasonably strong shock wave with boundary-layer 
separation near the trailing edge, the shock-wave location and the mean and 
rms pressure distributions were well predicted except over the last 20% of 
the airfoil. At a free-streara Mach number of 0.88, where the characteristic 
features of the steady flow are a shock wave sufficiently strong to induce a 
boundary-layer separation that extends from the base of the shock wave down- 
stream into the wake, the predictions were not as accurate as they were at 
Mach number 0.81, particularly for the chordwise location of the shock wave. 

At M.X, = 0.88 the predicted shock-wave location was 13% downstream of the 
experimental location and at M« “ 0.81 less than 2% downstream. This is 
believed to result frem a greater effect on the cemputations of the mismatch 
in the boundary conditions when the sonic line in the flow field about the 
airfoil extends closer to the wind-tunnel walls at the higher free-stream 
Mach number. 

Experimental unsteady flows were observed over a Mach-number range from 
0.82 to 0.86 and were documented in detail at M<» * 0.85. These flows are 
characterized by periodic oscillations of the shock-wave location and inten- 
sity and in the location of boundary- layer separation on opposite sides of the 
airfoil 180° out of phase. Computations at a free-stream Mach number of 0.83 
predicted the identical features. Locations of the shock waves and their 
intensit and frequency were accurately predicted. Time histories of the 
fluctuati ig pressures were in good agreement with experiment at all chordwise 
stations except for 0.5 < x/c < 0.7. The computed mean and rms pressure 
distributions and time histories of the unsteady lift and pitching moment were 
in good agreement with experiment. 

The results of comparisons of the predicted and experimental flow fields 
in the present study are sufficiently encouraging that in the near future a 
computational capability for predicting the Reynolds number, Mach number, and 
angle of attack at which airfoils can be expected to experience unsteady 
forces and for predicting the magnitude of those forces can be anticipated. 
Just how far into the future depends on progress in developing improved tur- 
bulence models for unsteady flows and on the availability of well- documented 
experiments — particularly for the boundary conditions. 
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